String breaking with Wilson loops? 



Slavo Kratochvila*** and Philippe de Forcrand^ 

'"Institute for Theoretical Physics, ETH Ziirich, CH-8093 Ziirich, Switzerland 
''Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

A convincing, uncontroversial observation of string breaking, when the static potential is extracted from Wilson 
loops only, is still missing. This failure can be understood if the overlap of the Wilson loop with the broken string 
is exponentially small. In that case, the broken string ground state will only be seen if the Wilson loop is long 
enough. Our preliminary results show string breaking in the context of the 3d SU{2) adjoint static potential, 
using the Liischer-Weisz exponential variance reduction approach. As a by-product, we measure the fundamental 
SU (2) static potential with improved accuracy and see clear deviations from Casimir scaling. 



1. INTRODUCTION 

The breaking of a long flux tube between two 
quarks into a matter- antimatter pair is one of the 
most fundamental phenomena in QCD. Three ap- 
proaches have been used to measure the static po- 
tential: (i) Correlation of Polyakov Loops, at fi- 
nite T [0 . (ii) Variational Ansatz using two types 
of operators: One for the string- like state and one 
for the broken string state . (iii) Wilson loops 

Hi- 

String breaking has been seen using the first 
two methods, but no clear signal using the third 
one has been observed. The failure of the Wil- 
son loop method seems to be mainly due to the 
poor overlap of the operators with the broken 
state. The strong coupling model |^ uses the 
heavy quark expansion to show that there is an 
exponentially small overlap. Based on a topo- 
logical argument, suggests that there may be 
no overlap at all. We show that there is a small 
overlap but not nearly as small as predicted. 

2. STATIC POTENTIAL 

The static potential helps to characterize con- 
fining forces. Assume that a static charge and 
a static anticharge are separated by a distance 
R < Rb, where Rb is the string breaking dis- 
tance. If R is increased to a value larger than Rf,, 
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a matter-antimatter pair will be created which 
screens the static charges. A further increase of 
R has no effect on the static potential - it remains 
« constant, this indicates string breaking. In our 
case, we deal with adjoint static charges which are 
screened by the gluons of the gauge field. The ob- 
ject adjoint-charge-gluon is called a gluelump. 

In a Hamiltonian formulation, a static charge 
and a static anticharge are created at Euclidean 
time 0, forming a state and annihilated at 
time T. This can be expressed as the Wilson loop 
W{R,T): 

{cl,\e-^"\cl>)=W{R,T) (1) 

On the other hand, one can expand the left side 
in (|l|) in the eigenbasis Vt^"^ of the Hamiltonian. 
At large T, only the ground state survives: 

^ |(M.(")|0)pe-™ (2) 

By using (Q) and (||), the static potential V{R), 
which is defined as the energy of the ground state 
Eq, can be connected with the Wilson loop: 

V{R) = - \im ^ log W{R,T) (3) 

This holds as long as the overlap ((/)|^'^*'^) does 
not vanish. 
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In section U, we fit the energy of tlie unbroken 
case using tlie Ansatz: 

V{R) ^Vo-^+aR (4) 
wliere a is tlie string tension. 



3. METHOD 

We are using a 3d-lattice with extent 48^ x 64 
at /3 = ^ = 6.0. The gauge group is SU{2). 
We measure adjoint static charges in order to 
avoid costly dynamical quarks. The price to pay 
is that the signal decays much faster because 
Vadj « l^/und (Casimir scaling). 

We make use of Liischer-Weisz exponential 
variance reduction This method general- 

izes Multihit [D from single links to link-link- 
correlators T{na). The possibility to use a hier- 
archical scheme allows us to measure long Wilson 
loops up to the desired accuracy. Using the no- 
tation introduced in we actually implement 
the following scheme: 



W{R,T)^-L{OUA- 



[T(na)][T((n+l)a)] 



[T((n+2)a)][T((n+3)a)] 



(5) 

The parameters are the following: The inner- 
most average [T(ria)] was calculated from 10 
sets of time-like links, each obtained after 10 
updates (1HB/40R). The second-level average 

[T(na)][T((n -I- l)a)] was calculated from 160 

averages of [T(na)] separated by 200 updates of 
the spatial links on the time-slice [{n + l)a]. Fi- 
nally, the third-level average was calculated from 
39 second-level averages separated by 200 updates 
of the spatial links on the time-slice [{n + 2)a]. 

To reduce contributions from excited states 
\]/("#o)^ we replace the simple spatial transporters 
L(0)c(7 ^C^Tps ■with "staples", which are 

constructed in the following way: (i) After each 
calculation of second level averages, I I , we form 

the smeared adjoint spatial links, , at the 

time-slice {n + 2)a. (ii) We multiply them with 



the second level averages to obtain the staples 

I ^ I K (iii) This procedure is repeated 39 times 

to provide an error reduction in the new spatial 
transporters. 

After processing 40 configurations as above, we 
obtain the following preliminary results. 

4. RESULTS 

We measure both the adjoint and the funda- 
mental static potential (Fig. 1). The spatial sep- 
aration R starts at 2 and is increased up to 12. In 
the adjoint case, the Ansatz (Q) works well for the 
energy V' (i?) of the unbroken string in the whole 
range of R. However, at i? > 10 the minimal 
energy remains constant, V{R > 10) = 1.95(10). 
Therefore, we observe string breaking at a dis- 
tance Rb « 10. 

The dotted fine at 2.06(1) represents twice the 
energy of a gluelump using a direct measurement 
Q on our lattice. This value is close to the 
groundstate potential for R > Rb (See Fig.l). 
Both Rb f 



and 2E{Qg) 




I, as the string breaking distance, 
|, as twice the gluelump-energy, 
are in agreement with the literature. 

In the fundamental case, the Ansatz (^) again 
works very well. An important issue is the ful- 
filment of the Casimir scaling law. At i? < 4, 
~ f holds as required by perturbation the- 
ory. However at large distances, the slopes of the 
potentials are in the ratio 2.31(1) < |. 
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Fig. 1. The adjoint and ^ fundamental static potentials 
V{R) vs R. The horizontal line at 2.06(1) represents twice 
the energy of a gluelump. 
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The value of the static potential V{R) at a 
fixed R, e.g. R — 12, is determined by the 
slope of the logarithm of the Wilson loop (^). In 
Fig. 2 it can be seen that a single exponential 
fit is not sufficient in the adjoint case. At small 
temporal extent T we get a larger slope than at 
large T. This can be explained as follows: At 
small T the signal is dominated by the unbroken 
string state. The broken string state can only 
be observed after T is large. The way to fit all 
data points is to use the double-exponential fit 
~ coe-^'(«)^ + cie-^-ElQs)?-. V'{R) corresponds 
to the energy of the unbroken string, E{Qg) to the 
energy of the gluelump. The ratio of the ampli- 
tudes is ^ 10'', the turning point is at Tt « 20. 
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Fig. 2. Adjoint Wilson loop W{R,T) vs T, for R = 12. 
The double-exponential fit matches all data points. Single- 
exponential fits (dotted lines) do not. The unbroken string 
state can be fitted by using data from T = 8 up to 20, the 
broken string state by using data from T = 20 up to 28. 
Note how small a signal can be measured. 

5. CONCLUSION 

In the context of the 3d SU{2) adjoint static 
potential, string breaking has been observed using 
Wilson loops only. 

It is a necessary condition that both R and T 
have to be larger than the string breaking dis- 
tance Rh. Note, however, that the variational ap- 
proach allows to bypass this requirement. 

The strong coupling model ^ is too pes- 
simistic. A necessary condition to be able to mea- 
sure the broken state is that the signal should 



be larger than the one from the unbroken state: 
^^^-2Mq,t ^ coe~^^^'>'^. The heavy quark ex- 
pansion results in the ratio ^ ~ ^MQgR^ which 
would imply T 3> 43 at i? = 12, whereas in our 
simulation the turning point was Tt w 20. 

A similar observation using Wilson loops in full 
QCD should be diflicuh but feasible 
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